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Light quark masses are computed for Kogut-Susskind fermions by evaluating non-perturbatively the renormal- 
ization factor for bilinear quark operators. Calculations are carried out in the quenched approximation at f3 = 6.0, 
6.2, and 6.4. For the average up and down quark mass we find mj^g-(2GeV) = 4.15(27)MeV in the continuum 
limit, which is significantly larger than 3.51(20)MeV (q* = 1/a) or 3.40(21)MeV (q* = 7r/a) obtained with the 
one-loop perturbative renormalization factor. 



1. Introduction 

Light quark masses are important unknown pa- 
rameters of the standard model, and a number of 
lattice QCD calculations have been carried out 
to evaluate quark masses employing the Wilson, 
clover or Kogut-Susskind (KS) fermion action [Q. 
Among the results, those with the KS action ap- 
pear more accurate than others because of small 
lattice discretization errors and small statistical 
errors. 

A worry with the KS result, however, has been 
that the employed one-loop renormalization fac- 
tor takes a large value of rj 2 in the range of 
(3 studied, calling into question the viability of 
perturbation theory. In this article we report a 
study to circumvent this problem: we calculate 
the renormalization factor of bilinear quark op- 
erators for the KS action non-perturbatively us- 
ing the method of Ref. M developed for the Wil- 
son/clover actions. This calculation is carried out 
in quenched QCD at j3 = 6.0, 6.2, and 6.4 on an 
32 4 lattice. The results, combined with our pre- 



vious calculation of bare quark masses Q], lead 
to a non-perturbative determination of the light 
quark mass. 

2. Method 

The renormalization factor of a bilinear opera- 
tor O is obtained from the amputated Green func- 
tion. 



T {p) = S(p)- l (0\<l>(x>)Ofo)\0}S(x>r 



(1) 



where the quark two-point function is defined by 
S(p) = (0\4>(p)4>(p)\0). The quark field (f>(p) with 
momentum p is defined from the original one- 
component field x( x ) by 4>a{p) — J2 y ex P( — W ' 
y)x(y + aA )i where y p = 2an M , p^ = 2ir/(aL)n fl 
(n M = [-L/A, L/4)) and A„ = [0, 1]. 

The renormalization condition imposed on 
To(p) is given by 



T {p) = Z^p)Zo{p)T% ] 



(2) 



where I^ is the amputated Green function at 
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tree level and Z^p) is the wave function renor- 
malization factor which can be calculated by the 
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Figure 1. The scalar renormalization factor 
Zs(p) and that for the pseudoscalar Zp(p). 

condition Zy (p) = 1 for the conserved vector cur- 
rent. 

The relation between the bare operator on the 
lattice and the renormalized operator in the con- 
tinuum takes the form, 

^ms(^) = U m (v,p)Z m (p)/Z (p)O lat -(a) (3) 

where U^m([i,p) is the renormalization-group 
running factor, and Zy^(p) is the matching fac- 
tor from the RI scheme defined by (§) to the MS 
scheme, calculated perturbatively in the contin- 
uum. For the light quark mass we apply relation 
([|) in the scalar channel in the chiral limit. 

We use a source in momentum eigenstate to 
evaluate quark propagators. This results in very 
small statistical errors of O(0.1%) in the Green 
functions. 

The external momentum p should be taken in 
the range Aqcd <C p <C 0(1/ a) in order to keep 
under control higher order effects in continuum 
perturbation theory, non-perturbative hadroniza- 
tion effect on the lattice, and discretization errors 
on the lattice. In this work we choose 15 mo- 
menta in the range 0.038553 < (ap) 2 < 1.9277 
for all values of (3. 

3. Result 

In Fig. [l]we compare the scalar renormalization 
factor Zs(p) with that for pseudoscalar Zp(p) for 
three values of bare quark mass am at = 6.0. 
From chiral symmetry of KS fcrmions, we expect 
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Figure 2. M(p) in the chiral limit. 

naively Zg(p) = Zp(p) for all momenta p in the 
chiral limit. Clearly this relation does not hold 
with our result toward small momenta, where 
Zp(p) rapidly increases as m — > 0, while Zg(p) 
does not show such a trend. 

To understand this result, we note that chiral 
symmetry of KS fermion leads to the following 
identities : 

Z s (p)-Z<f,(p)=dM(j>)/dm 

Z P (p) ■ Z+(p) = M(p)/m (4) 

with M(p) = Tr[5(p) -1 ]. In Fig. | M(p) in the 
chiral limit obtained by a linear extrapolation in 
m is plotted. It rapidly dumps for large momenta, 
but takes large values in the small momentum re- 
gion. Combined with (Q) this implies that Zp(p) 
diverges in the chiral limit for small momenta, 
which is consistent with the result in Fig. |]. 

The function M(p) is related to chiral conden- 
sate as follows : 

P 



^fC^(pf + M(pf 



(5) 



where C^p) = -iTr[( 1 ^I)S(p)- 1 }/ cos(p^a). A 
non- vanishing value of M(p) for small momenta 
would lead to a non-zero value of the condensate. 
Therefore the divergence of Zp(p) near the chiral 
limit is a manifestation of spontaneous symme- 
try breakdown of chiral symmetry; it is a non- 
perturbative effect arising from the presence of 
massless Goldstone boson. 

We expect this non-perturbative effect to affect 
the scalar renormalization factor Zg(p) much less, 
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Figure 3. The ratio m^(fx)/m at /i = 2GeV. 
For each (3 the filled data points are used for linear 
extrapolation in (ap) 2 . 

since the scalar operator can not interact directly 
with the pseudoscalar meson. Indeed the quark 
mass dependence is quite small as we have seen 
in Fig. 0. 

In Fig. |^ we show the momentum dependence 
of the ratio m m (fi)/m = U m {p,,p)Z m (p)Z s (p) 
calculated in the chiral limit where we set /i = 
2GeV and use the three-loop formula ||] for U-^r 
and Z-g^r. While the ratio should be indepen- 
dent of the quark momentum p, our results show 
a large momentum dependence which is almost 
linear in (ap) 2 for 0.6 < (ap) 2 . 

A natural origin of the linear dependence on 
(ap) 2 is the lattice discretization error of the 
scalar operator, which differs by terms of 0(a 2 ) 
from that of the continuum for the KS fermion. 
We then remove this error from the renormaliza- 
tion factor by a linear extrapolation in (ap) 2 to 
(ap) 2 = 0. In Fig. || the fitting lines are plotted, 
where filled data points are used for the linear ex- 
trapolation. For comparison, the ratio calculated 
with the one- loop value equals 1.867, 1.877, and 
1.871 for (3 = 6.0, 6.2 and 6.4 at q* = I /a. Hence 
one-loop perturbation theory underestimates the 
ratio by 40% to 20%. 

Our final results for the averaged up and down 
quark mass at /i = 2GeV are shown in Fig. |J by 
filled symbols. Here we use the JLQCD results 
for bare quark mass Q. The values are substan- 
tially larger than those obtained with one-loop 
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Figure 4. The final results of the light quark 
mass at \i = 2GeV. 

perturbation theory (open circles for q* = X/a 
and squares for q* — ir/a). Furthermore they 
exhibit a significant a 2 dependence, which we as- 
cribe to the discretization error of the quark mass 
itself. Making a linear extrapolation in a 2 , our fi- 
nal result in the continuum limit is given by 

m M g(2GeV) = 4.15(27)MeV. (6) 

This value is 20% larger than the perturbative 
estimates : 3.51(20)MeV for q* = 1/a and 
3.40(21)MeV for q* = vr/a. 
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